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Contrarily to general believe, a first-order cosmological perturbation theory based on Einstein's 
General Theory of Relativity explains the formation of massive primeval stars in a flat Friedmann- 
Lemaitre-Robertson- Walker universe after decoupling of matter and radiation, whether or not Cold 
Dark Matter is present. The growth rate of a density perturbation depends on the heat loss of a 
perturbation during the contraction, but is independent of the particle mass. The relativistic Jeans 
mass does depend on the particle mass. If the Cold Dark Matter particle mass is equal to the proton 
mass, then the relativistic Jeans mass is equal to 3500 solar masses, whereas the classical Jeans mass 
is a factor 145 larger. 
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I. INTRODUCTION 

A manifestly covariant gauge-invariant cosmological 
perturbation theory for flrw universes based on the 
Theory of General Relativity combined with Thermody- 
namics and a realistic equation of state for the pressure 
p = p(n, e) has been developed in a foregoing article [1]. 
In fact, in this article the pioneering work of Lifshitz and 
Khalatnikov [2, 3] has been redone, but now based on 
two newly introduced gauge invariant quantities. These 
quantities, which we baptized sf\ and nfl^ turned out 
to be the energy density and the particle number den- 
sity perturbations. Indeed, taking the non-relativistic 
limit v/c — > the complete set of relativistic perturba- 
tion equations reduce to the complete set of Newtonian 
equations 

WW = 4rf^, efi ) (x) = nf 1) (x)mc?, (1) 

thus identifying ef[ } and the real energy density 

and particle number density perturbations. 

In the present article it will be demonstrated that 
the formation of primordial stars, the so-called (hypo- 
thetical) population in stars [4, 5], can be understood 
with the revised Lifshitz-Khalatnikov perturbation the- 
ory, whether or not Cold Dark Matter (cdm) is present. 

A. Former Results: Standard Perturbation Theory 

It is generally accepted that in a universe filled with 
only 'ordinary matter,' i.e., elementary particles and pho- 
tons but no CDM, linear perturbation theory predicts a 
too small growth rate to explain the formation of struc- 
ture in the universe. The reason put forward in the liter- 
ature on structure formation is that in the linear phase 
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of the growth of an adiabatic relative density perturba- 
tion S(t, x) in the era after decoupling of radiation and 
matter, given by 

/ t \ 2/3 

S(t) = 6(t dec ) — , t dcc <t<t p , (2) 

\ tdcc / 

is insufficient for relative density perturbations as small 
as the observed initial value S(tdec) ~ 10 -5 to reach the 
non-linear phase for times t <t p , where i p = 13.75 Gyr, 
the present age of the universe, and tdcc = 381 kyr, the 
time of decoupling of matter and radiation [6, 7]. This 
generally accepted conclusion follows from the standard 
evolution equation for density perturbations in a uni- 
verse which is after decoupling of matter and radiation 
assumed to be filled with a perfect and pressure-less fluid 
(usually referred to as 'dust') with equations of state for 
the energy density e and pressure p 

e = nmc 2 , p = 0, (3) 

where n is the particle number density and m the particle 
mass. 

Before decoupling, Thomson scattering between pho- 
tons and electrons and Coulomb interactions between 
electrons and baryons were so rapid that the photons and 
baryons are tightly coupled so that the photon-baryon 
system behaves as a single fluid. The standard pertur- 
bation theory predicts that density perturbations in this 
baryon-photon fluid oscillates with a constant amplitude 
and thus do not grow at all before decoupling. Since CDM 
is supposed to be electrically neutral, it is not linked by 
Coulomb interactions to the baryon-photon fluid. There- 
fore, researchers [8, 9] in the field of structure formation 
have assumed in their simulations that CDM would have 
already clustered before decoupling and thus would have 
formed seeds for baryon contraction after decoupling. If 
this would be true, then the slow growth (2) would be 
sufficient to explain structure in the universe. Thus, the 
mechanism of structure formation relies heavily on the 
particular property of CDM before decoupling, namely 
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that density perturbations in CDM are not electrically 
coupled to perturbations in the total energy density. 



B. Results from the Revised Lifshitz-Khalatnikov 
Perturbation Theory 

However, in a foregoing article [1] it has been demon- 
strated that in the radiation-dominated era density per- 
turbations in ordinary matter and CDM are both gravi- 
tationally coupled to density perturbations in the total 
energy density. Moreover, it has been found that small- 
scale density perturbations oscillate in the radiation- 
dominated era with an increasing amplitude, propor- 
tional to i 1 / 2 . These energy (i.e., radiation, ordinary 
matter and CDM tightly coupled together) density per- 
turbations will form the seeds for star formation after 
decoupling. They manifest themselves as small temper- 
ature fluctuations (24f) in the cosmic background radia- 
tion. 

Finally, it has been shown that neglecting the kinetic 
energy density ^nk^T with respect to the rest energy 
density nmc 2 yields in the perturbed universe the non- 
relativistic limit with gfj, = nfl } mc 2 and p {0) = p%. = 0, 
(3), implying that 5 £ = S n . Perturbations described by 
equations of state (3) are adiabatic. The growth is in this 
case given by the standard expression (2), which is far 
too slow to account for structure in the universe. Since 
adiabatic density perturbations cannot lose their internal 
energy to their environment, they grow only under the 
influence of gravitation. This explains their slow growth. 

In the dark ages of the universe (i.e., the epoch be- 
tween decoupling and the ignition of the first stars) a 
density perturbation from which stars will eventually be 
formed should have initially a somewhat smaller internal 
pressure than its environment and has to lose some of 
its heat energy in order to grow faster than given by the 
standard growth rate (2). It has been established [1] that 
in a non-static universe density perturbations described 
by an equation of state p — p(n, e) are diabatic, whereas 
perturbations described by (3) are adiabatic. Therefore, 
the evolution of density perturbations should be studied 
by using a realistic equation of state for the pressure of 
the form p = p(n, e), so that next to gravitational forces, 
also the heat exchange of a perturbation can be incorpo- 
rated into the perturbation theory. In fact, incorporating 
the realistic equations of state (4) , yields in the final dy- 
namical perturbation equation (8a) a source term and 
an evolution equation (8b) for this source term. Using 
the combined First and Second Laws of Thermodynam- 
ics dE = TdS — pdV + jtxdJV the source term of equation 
(8a) can be identified with the entropy of a perturbation. 
From equation (13) one may infer that the pressure term 
and the entropy term are of the same order of magni- 
tude. This yields in the early stages of the contraction 
of a perturbation a somewhat larger growth rate than in 
the adiabatic case (2). This faster growth is just enough 
for density perturbations with initial values as small as 



S n ~ S e < 10 5 to reach the non-linear regime within 
10 2 -10 3 Myr. 

II. OUTLINE 

Since CDM and ordinary matter particles behave grav- 
itationally in exactly the same way and since the mass 
of a CDM particle is as yet unknown, we assume that the 
CDM particle mass is approximately equal to the proton 
mass. 

After decoupling the cosmic fluid is a mixture of 
baryons (protons) and CDM. This mixture can be con- 
sidered as a non-relativistic monatomic perfect gas with 
equations of state for the energy density and the pressure 

e(n,T) — nmc 2 + ^nk^T, p(n,T) — nk^T, (4) 

where fee is Boltzmann's constant, m the mean particle 
mass, and T the temperature of the matter. It is assumed 
that m = tojj = ?7icDM> where mu is the proton mass 
and mcDM the mass of a CDM particle, implying that 
mc 2 3> k^,T throughout the matter-dominated era af- 
ter decoupling. Therefore, one may neglect the pressure 
nk#T and kinetic energy density ^nk^T with respect to 
the rest-mass energy density nmc 2 in the unperturbed 
universe: the kinetic energy of the particles in the uni- 
verse has a negligible influence on the global evolution 
of the universe. Thus, as is well-known, the global prop- 
erties of the universe after decoupling are very well de- 
scribed by a perfect and pressure-less fluid ('dust'), de- 
scribed by equations of state (3). 

At the moment of decoupling of matter and radiation 
photons could not ionize matter any more and the two 
constituents fell out of thermal equilibrium. As a con- 
sequence, the pressure drops from a very high radiation 
pressure p = T^aB,T^ just before decoupling to a very low 
gas pressure p = nk^T after decoupling. This fast and 
chaotic transition from a high pressure epoch to a very 
low pressure era may result locally in large relative pres- 
sure perturbations. These pressure perturbations will 
be taken into account by incorporating the equations of 
state (4) and their perturbed counterparts 

S n - 5 £ re -^. kBT( -° } s T , 5 p = 5 n + S T , (5) 
2 mc z 

into the new perturbation theory. In the expressions (5), 
T (0) is the background matter temperature and S n , S e , 
St and S p are the relative perturbations in the parti- 
cle number density, the energy density, the matter tem- 
perature and the pressure, respectively. The influence 
of pressure perturbations on the growth of small den- 
sity perturbations can only be investigated by using the 
revised Lifshitz-Khalatnikov perturbation theory devel- 
oped in [1], since this theory not only has an evolution 
equation (8a) for S £ , but (in contrast to all former per- 
turbation theories) also an evolution equation (8b) for 
the difference 8 n — S £ . This proves to be crucial for the 
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understanding of star formation in the early universe. 
Although in a linear perturbation theory \S P \ < 1 and 
\$t\ < 1, the initial values of these quantities are, accord- 
ing to (5), not constrained to be as small as the initial 
values 



5e(td ec ,q) « S n (t dec ,q) < 10" 



(6) 



as is demanded by WMAP-observations [6, 7]. Since the 
gas pressure p — nk^T is very low, its relative perturba- 
tion S p = p^)/p( 0) and, accordingly, the matter tempera- 
ture perturbation St = T^ } /T (Q) could be large. 



III. RESULTS 

Just after decoupling, ordinary matter is mixed with 
CDM. It is found that the growth rate is independent of 
the particle mass, i.e., the gravitational mechanism for 
star formation works equally well with or without CDM. 

It will be shown that just after decoupling at z = 1091 
negative relative matter temperature perturbations as 
small as —0.5% yields massive stars within 13.75 Gyr. 
The very first stars, the so-called Population in stars, 
come into existence between 10 2 Myr and 10 3 Myr. The 
star masses are in the range from 4 x 10 2 M© to 10 5 M Q , 
with a peak around 3.5 x 10 3 M©. Density perturbations 
with masses smaller than 3.5 x 10 3 M© become non-linear 
at later times, because their internal gravity is weaker. 
On the other hand, density perturbations with masses 
larger than 3.5 x 10 3 M© enter the non-linear regime also 
later, since they do not cool down so fast due to their 
large scale. The mass 3.5 x 10 3 M© corresponds initially 
to a scale of 6.2 pc. These conclusions are outlined in 
Figure 1. From this figure it follows that the growth 
rate rapidly decreases for perturbations with masses be- 
low 3.5 x 10 3 M©. Therefore, the peak values in Figure 1 
can be considered as the relativistic counterparts of the 
classical Jeans mass. However, the Jeans mass does de- 
pend on the particle mass: heavier particles yield lighter 
primordial stars. 



IV. BASIC EQUATIONS 

For the equations of state (4) the background equa- 
tions for a flat (k = 0) flrw universe with a vanishing 
cosmological constant (A = 0) reduce to 



3H 2 



-me, 



-3Hn, 



(7) 



where it is used that w = p (a )/e {0) < 1, so that the 
background pressure p m can be neglected with respect to 
the background energy density e (0) . An overdot denotes 
differentiation with respect to ct, and k = 8irG/c . 

It has been shown in a foregoing article [1] that for 
equations of state (4) and their perturbed counterparts 



(5) the perturbation equations reduce to 
,V 2 



S £ + 3H5 
1 d 



2V %A 
'3^ (<5 " 



cdt 



(S n - S € ) = -2H (S n - S E ) . 



(8a) 
(8b) 



The quantity /3(t) defined by /3 = v^mT^w i s i to a g°°d 
approximation, given by 



v s (t) 



5 k B T m (t) 



T m cx a 2 , (9) 



with v s the adiabatic speed of sound and T (0) the matter 
temperature. Equation (8b) implies 



S n — S e oc a 2 , 



(10) 



where it is used that H = a/a, with a(t) the scale factor 
of the universe. Combining (9) and (10) one gets from (5) 



5 T (t,x) « 5<r(*o,a0, 



(11) 



to a very good approximation. Using the well-known 
solutions of the background equations (7) 



H cx t , £ (0) (xt 2 , n (0) cx t 2 , a oc t 



2/3 



(12) 



and substituting 8(t,x) = S(t, q) exp(iq ■ x), equations 
(8) can be combined into one equation 



5" + -S' + 
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9 r 8 / 3 9r 2 



±]4n_ 

15 r 8 / 3 



S T (t ,q), (13) 



where r = t/tg and a prime denotes differentiation with 
respect to r. The parameter ^ m is given by 



2ir 1 v s (t ) 
A H(t ) c 



, A = Xa(t ), 



(14) 



where Ao = 27r/|q | is the scale of a perturbation at 
time £o- The constant fi m can be expressed in observ- 
able quantities. To that end we use that the redshift z(t) 
as a function of the scale factor a(t) is given by 



z(t) = 



o(t) 



(15) 



where a(t p ) is the present value of the scale factor. For a 
flat flrw universe one may take a(t p ) — 1. Using (12), 
(15) and T (0) cx a -2 , we get 



2tt 



Mm = 



5 fc B r ( o)(tdcc) 

3 m 



A H(t p )[z(t dec ) + 1] ^z(t a ) + r 



(16) 



where to is the time when a perturbation starts to con- 
tract. This expression is invariant under the replacement 
m — > am and Aq — > Xo/y/a, for some constant a > 0. 
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This implies that a perturbation S E with initial scale Ao 
in a cosmic fluid with mean particle mass m, evolves 
in exactly the same way as a perturbation with initial 
scale Xo/t/cx in a fluid with mean particle mass am. In 
other words, the growth rate is independent of the par- 
ticle mass. 

The mass M(i ) of a spherical density perturbation 
with radius ^Aq is given by 



M(t ) = y (iA ) 3 n (0) (£ )m. 



(17) 



The particle number density n (0) (io) can be calculated 
from its value nf 0) (t cq ) at the end of the radiation- 
dominated era. By definition, at the end of the radiation- 
domination era the matter energy density n {0) mc 2 equals 
the energy density of the radiation: 

n (0) (t oq )mc 2 = a B T ( o )7 (ioq)- (18) 

Since n (0) oc a -3 and T (0)7 oc a -1 , we find, using (15), the 
particle number density at time to 



a BT (0) Jt p 



mc* 



[z(t cq ) + 1] [z(t ) + 1] 



(19) 



Combining (17) and (19), we get for the mass of a spher- 
ical density perturbation 



M(t 



47T 



.3 flB^o^^p) 



[z(i eq ) + l] [>(to) + l] 



(20) 

With the help of this expression the initial scale Ao of a 
perturbation is related to its mass at the initial time to- 
The influence of the mean particle mass m on the mass 
M(to) of a primordial star can be studied by replacing 
m by am and Ao by Ao/v^ {a > 0) in expressions (16)- 
(19). It is found from (20) 



M(t ) cx a~ 3/2 . 



(21) 



In other words, the heavier the particles in the universe, 
the lighter the primordial stars. For example, if a = 10 
then the mean particle is m = 10?tih, implying that the 
Jeans mass in Figure 1 becomes approximately 1O 2 M . 

Finally, the influence of the initial time on a star mass 
is determined. It follows from (16) that 



A oc [zfa) + 1] 



-1/2 



implying with (20) that 

M (to) cx [z(t ) + if 2 



(22) 



(23) 



Thus, the later a perturbation starts to contract, the 
smaller the mass will be. For example, if a perturba- 
tion starts to contract at z(t a ) = 1, then the Jeans mass 
in Figure 1 will be Mj(t ) ~ 0.27 M . 



V. INITIAL VALUES FROM WMAP 

The physical quantities measured by wmap [6, 7] and 
needed in the present theory of primordial star formation 
are the redshifts at matter-radiation equality and decou- 
pling, the present values of the Hubble function and the 
background radiation temperature, the age of the uni- 
verse and the fluctuations in the background radiation 
temperature: 

z(t cq ) = 3196, (24a) 
z(t dcc ) = 1091, (24b) 
off(tp) = H{t p ) = 71.0km/scc/Mpc, (24c) 
r (0)7 (t p ) = 2.725 K, (24d) 
t p = 13.75 Gyr, (24e) 

M*dec) <nr 5 . (24f) 

At decoupling the matter temperature is equal to the 
radiation temperature. The latter can be calculated from 
the fact that T (0)7 cx a" 1 and the quantities (24b) and 
(24d). Using (15), one finds for the matter temperature 
at decoupling 



T'(o)(tdoc) — 7 1 (0)7 (tdoc) = 2976 K. 



(25) 



Substituting the observed values (24b) and (24c) into 
(16), one gets, using also (25), 



Mb 



518.5 



AoaA(*o) + 1 : 



Aq in pc, 



(26) 



where we have used that 1 pc = 3.0857 x 10 16 m (1 pc = 
3.2616 ly). 

Finally, using that one solar mass is 1.98892 x 10 30 kg, 
we find from (24) that 



M(t ) = 1.148 x 10" 8 A 3 [z(t ) + 1] M . 



(27) 



The expression (27) will be used to convert the scale Ao 
(expressed in units of 1 pc) of a perturbation, which starts 
to contract at a redshift of z(to)i into its mass M(to) 
(expressed in units of the solar mass). 



VI. POPULATION III STAR FORMATION 

In this section the evolution equation (13) is solved 
numerically. To that end the differential equation solver 
lsodar with root finding capabilities is used. This solver 
is included in the package deSolve [10], which, in turn, 
is included in R, a system for statistical computation and 
graphics [11]. Star formation which starts at cosmologi- 
cal redshift z = 1091, i.e., at to = tdec, is investigated. 

The WMAP observations of the fluctuations in the back- 
ground radiation temperature yield for the fluctuations 
in the energy density and particle number density (6). In 
addition, it is assumed that 



4 (tdec, q) « 0, 



(28) 



5 



Star Formation starting at z = 1091 




j 

le+02 le+03 le+04 le+05 
Star Mass in Units of the Solar Mass 



le+06 



Figure 1. The curves give the redshift at which a linear per- 
turbation in the particle number density with initial values 
S„(tdec,q) ~ S e (t dcc ,q) w 10" 5 and 5' n (tdec,q) = start- 
ing to grow at an initial redshift of z(ta cc ) = 1091 becomes 
non-linear, i.e., {, « { e w 1, During the evolution we have 
S p (t,q) — SritdecQ) + S n (t,q). The numbers at each of the 
curves are the initial relative perturbations in the matter tem- 
perature 5r(tdec, q). For each curve, the Jeans mass (i.e., the 
peak value) is at 3.5 x 10 3 M Q . 



i.e., during the transition from the radiation-dominated 
era to the era after decoupling energy density perturba- 
tions are approximately constant with respect to time. 

Figure 1 has been constructed as follows. For each 
choice of ^r^dec, q) equation (13) is integrated for a large 
number of values for the initial scales Ao = Adcc, using 
the initial values (6) and (28). The integration starts at 
t = t/td cc = 1, i.e., at z = z(t dcc ) an d will be halted if 
either z = (i.e., r = [z(t dec ) + l] 3 / 2 ), or S e (t,q) = 1 
has been reached. One integration run yields one point 
on the curve for a particular choice of the scale Ad ec if 
S e (t, q) — 1 has been reached for z > 0. If the integration 
halts at z = and still 5 e (t, q) < 1, then the perturbation 
belonging to that particular scale Adcc has not yet reached 
its non-linear phase today, i.e., at t p = 13.75 Gyr. On the 
other hand, if the integration is stopped at 6 e (t,q) = 1 
and z > 0, then the perturbation has become non-linear 
within 13.75 Gyr. 

The above described procedure is repeated for 
S T (tdcc,q) in the range -0.005, -0.01, -0.02, -0.1. 
During the evolution, the relative pressure perturbation 
evolves according to (5) and (11): 

S p (t,q) = S T (t dec> q) + 5 n (t,q). (29) 

The fastest growth is seen for perturbations with a mass 



of approximately 3.5 x 10 3 M Q . This value is nearly in- 
dependent of the initial value of the matter tempera- 
ture perturbation Jr^deo <?)■ Even density perturbations 
with an initial relative matter temperature perturbation 
as small as <5r(^dec, q) = —0.5% reach their non-linear 
phase at z — 0.13 (T (0)7 = 3.1 K, t = 11.5 Gyr) pro- 
vided that its mass is around 3.5 x 1O 3 M . Perturba- 
tions with masses smaller than 3.5 x 10 3 M Q reach their 
non-linear phase at a later time, because their internal 
gravity is weaker. On the other hand, perturbations with 
masses larger than 3.5 x 10 3 M Q cool down slower be- 
cause of their large scales, resulting also in a smaller 
growth rate. Since the growth rate decreases rapidly 
for perturbations with masses below 3.5 x 1O 3 M , the 
latter mass will be considered as the relativistic coun- 
terpart of the classical Jeans mass. This mass corre- 
sponds to a Jeans scale of 6.2 pc « 20 ly. This scale is 
much smaller than the horizon size at decoupling, given 
by du(tdec) = 3ci d cc ~ 350 kpc. 



VII. HEAT LOSS DURING CONTRACTION 

In this section the heat loss of a density perturbation 
during its contraction is calculated. To that end the com- 
bined first second law of thermodynamics (24) in Ref. [1] 
is rewritten in the form 



1 (0)*(1) 



-(Sn-5 e ) 



Pto) 



(30) 



where it is used that w = p {0) /e {0) . Substituting expres- 
sions (4) and (5) into (30) and using also (11), one finds 
the entropy per particle of a density perturbation: 

s&(t,a;)« tk B [3S T (to,x)-25 n {t,x)], (31) 

where it is used that mc 2 >• k B T (0) . For all values of 
tir^o!^) in Figure 1 and initial values (6) the entropy 
perturbation is negative, sfl } < 0. Since for growing per- 
turbations one has 8 n > the entropy perturbation de- 
creases, i.e., sfn = — k B 8 n < 0, during contraction. This 
implies that a growing perturbation loses a part of its in- 
ternal energy to its environment. This is to be expected, 
since a local density perturbation is not isolated from its 
environment. Only for an isolated system the entropy 
never decreases. 



VIII. CLASSICAL JEANS MASS 

In this section the classical Jeans mass, derived from 
the Newtonian theory of gravity, is compared with the 
relativistic Jeans mass which follows from the revised 
Lifshitz-Khalatnikov perturbation theory. 

The classical Jeans length at time £q is given by [12] 



Aj(to) = ^s(£o) 



Gn (0) (t )m' 



(32) 
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Using (9), (17) and (19) one finds for the classical Jeans 
mass at decoupling: 

MjM«5.1x 1O 5 M , (33) 

where the wmap values (24) have been used. The clas- 
sical Jeans mass (33) is much larger than the relativis- 
tic Jeans mass, 3.5 x 10 3 M which follows from the re- 
vised Lifshitz-Khalatnikov perturbation theory. This dif- 
ference is due to the fact that in the classical perturba- 
tion theory based on the equations of state (3) the heat 
loss of a perturbation is not taken into account, whereas 
the effect of heat loss on the growth of a perturbation 
is included in the revised Lifshitz-Khalatnikov perturba- 
tion theory based on the equations of state (4). In other 
words, since a perturbation loses some of its energy, grav- 
ity can be somewhat weaker to make a perturbation con- 
tract. The classical Jeans mass (33) corresponds to a 
classical Jeans scale of 32.4 pc s» 106 ly. Just as is the 
case for yU m (16), the expression (32) is invariant under 
the replacement m — > am and Aj — > Xj/y/a, for some 
constant a > 0. As a consequence, the classical Jeans 
mass Mj is proportional to a -3 / 2 , just as in the rela- 
tivistic case (21). 

Finally, the classical Jeans mass of a perturbation 
starting at z(to) = 1 follows from (23) and (33). It is 



found that Mj(i ) « 40 M Q . 



IX. CONCLUSION 

Three important conclusions can now be drawn. 
Firstly, there is no need to make use of alternative grav- 
itational theories: the Theory of General Relativity ex- 
plains the formation of massive primordial stars in our 
universe. In other words, Einstein's gravitational theory 
not only describes the global characteristics of the uni- 
verse, but is also locally successful. Secondly, although 
there is strong evidence for the existence of CDM [13], it is 
not needed for the formation of primeval stars. Finally, 
is has been demonstrated that not only for large-scale 
perturbations one should use the theory of relativity, but 
also for small-scale perturbations: because of the spu- 
rious gauge modes present in the Newtonian theory of 
gravity, it fails to predict primordial stars. 
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